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The testing of general relativity at cosmological scales has become a possible and timely endeavor 
that is not only motivated by the pressing question of cosmic acceleration but also by the proposals of 
some extensions to general relativity that would manifest themselves at large scales of distance. We 
analyze here correlations between modified gravity growth parameters and some core cosmological 
parameters using the latest cosmological data sets including the refined Cosmic Evolution Survey 
3D weak lensing. We provide the parametrized modified growth equations and their evolution. 
We implement known functional and binning approaches, and propose a new hybrid approach to 
' evolve the modified gravity parameters in redshift (time) and scale. The hybrid parametrization 

/"*\ | combines a binned redshift dependence and a smooth evolution in scale avoiding a jump in the 

matter power spectrum. The formalism developed to test the consistency of current and future data 
CO ' with general relativity is implemented in a package that we make publicly available and call ISiTGR 

(Integrated Software in Testing General Relativity), an integrated set of modified modules for the 
publicly available packages CosmoMC and CAMB, including a modified version of the integrated Sachs- 
Wolfe-galaxy cross correlation module of Ho et al and a new weak-lensing likelihood module for the 
' refined Hubble Space Telescope Cosmic Evolution Survey weak gravitational lensing tomography 

, data. We obtain parameter constraints and correlation coefficients finding that modified gravity 

parameters are significantly correlated with as and mildly correlated with Q m , for all evolution 
methods. The degeneracies between as and modified gravity parameters are found to be substantial 
for the functional form and also for some specific bins in the hybrid and binned methods indicating 
that these degeneracies will need to be taken into consideration when using future high precision 

?H ' data. 
> 

j^J PACS numbers: 95.36.+x,98.80.Es,98.62.Sb 
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I. INTRODUCTION 



CO 

> 

CO 

There has been a growing interest in testing general relativity (GR) at cosmological scales. This has been triggered 
by the quest to understand the origin of cosmic acceleration which is one of the biggest puzzles in cosmology today. 
Specifically, one would like to test whether cosmic acceleration is due to some unknown, " dark energy" component in 
0^ ' the Universe or rather an extension or modification to gravity physics on cosmological scales. The growth rate of large 
scale structure proved to be a solid discriminator between these two possible explanations for cosmic acceleration and 
. has been studied within this context by a number of papers. For example, some early papers explored inconsistencies 
between parameter constraints derived using observations that depended primarily on the growth rate of structure 
^ . versus constraints derived from observations that probed only the expansion history. The presence of such inconsis- 
tencies would call into question the underlying gravity theory (see [iHj]). These were followed by a number of papers 
using certain growth parameters that take distinct values for different gravity theories (see [5l445l|). In this second 
approach, some parameters characterizing the growth rate are introduced in order to parametrize deviations from 
GR. These parameters have known values in GR but different values in modified gravity theories. So one possible 
goal is to constrain these parameters, in addition to the usual cosmological parameters, and see if they are consistent 
or not with the GR values. Also, recently Basilakos et al. proposed using the evolution of the linear bias as a way to 
test for deviations from GR 46]. 

These approaches have been applied using current available data sets and future simulated ones and in order to 
constrain growth parameters related to modifications of the perturbed Einstein equations in an effort to look for any 
indications of deviations from GR, but so far, none have been detected |47l - l57j . 

In this paper, we study correlations between modified gravity parameters and core cosmological parameters. We 
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implement the evolution in time and scale of the modified gravity parameters using a functional form method and a 
binning method. We also propose and implement a new hybrid method that combines both. In the hybrid method, 
the evolution in time is represented in two redshift bins while the scale evolution follows a monotonic functional form. 
This provides a smooth evolution in scale combined with a binned redshift (time) dependence that was shown to 
be more robust than time functional forms as we discuss in the next section, see also [53|, [56J, l58| . We also describe 
the numerical framework that we introduce here as the package: Integrated Software in Testing General Relativity, 
ISiTGR (pronounced Is it GR and available publicly at |http : //www . utdallas . edu/~j d ossett/ isit g r/) . ISiTGR is 
an integrated set of modified modules for the publicly available packages CosmoMC 59] and CAMB [60]. It combines 
all the modifications to those packages and a modified version to the of the Integrated Sachs-Wolfe (ISW)-galaxy 
cross correlations module by Ho et al. [6l|, [62j to test GR. We also include our weak-lensing likelihood module for 
the recently refined Hubble Space Telescope (HST) Cosmic Evolution Survey (COSMOS) weak-lensing tomography 
analysis in [63[ which has also been modified to test GR, and a new baryon acoustic oscillation (BAO) likelihood 
module for the recently released WiggleZ Dark Energy Survey BAO measurement data [64| . 



II. PARAMETRIZING DEVIATIONS OF THE GROWTH EQUATIONS FROM GENERAL 

RELATIVITY 



A. Growth Equations in General Relativity 

In the conformal Newtonian gauge the perturbed Friedmann-Lemaitre- Robertson- Walker (FLRW) metric is written 



as 



ds z = a(rr[-(l + 2i))d-r + (1 - 2<j))dx l dx t ], (1) 

where a(r) is the scale factor normalized to 1 today, the cc^'s are the comoving coordinates, and r is conformal time. 
</> and ip are scalar potentials describing the scalar mode of the metric perturbations. 

Using the first-order perturbed Einstein equations, while working in Fourier k space, we can get two very useful 
equations that describe the evolution of the scalar potentials. The combination of the time-time and time-space 
equations gives the Poisson equation describing the potential 4>. Then, to relate the two potentials to one another we 
take the traceless, space-space component of these equations. Explicitly, these equations are 

k 2 (j> = -AirGa 2 ^ piAi (2) 

i 

fc 2 (V>-0) = -YlTtGat^p^l + w^au (3) 

i 

where pi and Cj are the density and the shear stress, respectively, for matter species, i. Aj is the gauge-invariant, 
rest-frame overdensity for matter species, i, the evolution of which describes the growth of inhomogeneities. It is 
defined by 

A 1= <5 1 + 3h|, (4) 

where % = a/a is the Hubble factor in conformal time, and for species i, Si = 5pi/p is the fractional overdensity and qi 
is the heat flux and is related to the divergence of the peculiar velocity, 9{, by 9{ = i + q ^. ■ Enforcing the conservation 
of energy momentum on the perturbed matter fluid, these quantities for uncoupled fluid species or the mass-averaged 
quantities for all the fluids evolve as described in 65] : 

SP 

S = -kq + 3{l + w)(/) + 'SH(w- — )S (5) 

Sp 

| = -W(l-3«;)? + ^5 + (l + «;)(^-ff). (6) 
k k op 

Above, w = p/p is the equation of state. Combining these two equations, we can express the evolution of A by 

A = 3(l + w)(<j> + Hi?j +3HwA~ k 2 + 3 (% 2 - H) | - 3H(l + w)a. (7) 

Equations ^ , Q , , and © are coupled to one another; combining them, along with the evolution equations for 
a(r), we can describe the growth history of structures in the Universe. 
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B. Modified Gravity Growth Parameters 

Parametrizing both modifications to Poisson's equation, ([2]), as well as the ratio between the two metric potentials 
<j> and if) in the perturbed FLRW metric (called gravitational slip by Caldwell et al. [Tlj ) has recently been the subject 
of a lot of the work on testing general relativity; see, for example, [TTI. I49i - f53| . The parameters we use in this paper 
to describe modifications to the growth [modified gravity MG parameters] are based upon those used in [Hoj . 

The parametrized modifications to the growth equations proposed by [5(| directly modify Eqs. ^ and ^ and 
make no assumptions as to the time when a deviation from GR is allowed. These modifications are as follows: 

k 2 (f> = -4irGa 2 J2p t A t Q (8) 

i 

fc 2 (V>-i?0) = -127rGa a £^ pj(l + «;<)£>•< Q, (9) 

i 

where Q and R are the MG parameters. The parameter Q represents a modification to the Poisson equation, while 
the parameter R quantifies the gravitational slip (at late times, when anisotropic stress is negligible, R = if) /4>). In our 
code, rather than using the parameter R which is degenerate with Q, we instead use the parameter V = Q(l + R)/2 
as suggested in (this parameter is equivalent to the parameter S in [IH [53[ or Q of [Hj]). Combining Eqs. ® 
and ©, we arrive at the second modified growth equation used in this paper: 

k 2 {if) + (f>) = -SirGa 2 ^ p 4 A, V - YlnGo 2 ^ p 4 (l + wjeri Q. (10) 

i i 

So, the modified growth equations are (JSJ and (|10p. and Q and V are now the MG parameters. As discussed in our 
previous work [561 ] . this approach of using the parameter T> instead of R is also useful because observations of the 
weak-lensing and ISW are sensitive to the sum of the metric potentials <f> + if) and its time derivative respectively. 
Thus observations are able to give us direct measurements of this parameter. 



C. Different Approaches to Evolving Modified Growth Parameters in Time and Scale 

To date, there have primarily been two approaches to evolving the MG parameters in time and scale; one using 
a continuous functional form and the other based on binning. We implement in ISiTGR the two approaches and, 
additionally, a new hybrid approach, as we explain below. 

The first approach involves defining a functional form for each parameter that allows it to evolve monotonically 
in both time and scale. This allows one to make no assumptions as to when a deviation from general relativity is 
allowed. Such an approach was taken in for example 501. In that work the functional form, 

X(k,a)= [x e~ fe/fcc +X 00 (l-e- k/k ")-l\a s + 1, (11) 

was assumed, where X denotes either Q or R in Eqs. © and ©■ Thus a total of six model parameters are used to 
test GR: Qq, Rq, Qoo, Roc, k c , and s. The parameters s and k c parametrize time and scale dependence respectively, 
with GR values s = and k c — oo. Qo and i?o are the present-day superhorizon values while Qoq and R^ are the 
present-day subhorizon values of the Q(k,a) and R(k,a) , all taking GR values of 1. 

In the second approach, instead of evolving each of the parameters assuming some functional form, one can bin the 
MG parameters. This approach allows the parameters to take on different values in predefined redshift and scale bins. 
This technique was used in, for example, (52l . [531 ] . In those works two redshift and two scale bins were defined and for 
redshifts above a certain critical redshift, GR was assumed to be valid. In each bin the parameters were allowed to 
take on different values resulting in a total of eight model parameters used to test GR. 

The third approach that we propose here is a hybrid one where the evolution in redshift (or time) is binned into 
two redshift bins, but the evolution in scale evolves monotonically in the same way as the functional form above. Our 
motivation for this hybrid binning approach is that it takes advantage of an evolution in scale that is not so abrupt 
as that in the traditional binning method, while still taking advantage of a redshift (time) dependence expressed in 
the form of bins, which was shown to be more robust than time functional forms (53l. [5a.[58j . For example, in [56j | . we 
found that binning methods do not display the extent of tensions between the MG parameters (preferred by different 
data sets) as in the functional form method, where tensions are exacerbated by the chosen functional form. Also, in 
Ref. 58], the author found similar to what we noticed and that the constraints on MG parameters depend strongly 
on the parameter s (the scale factor exponent in the functional parametrization) . They further looked at ways to 
remove this strong dependence on the parameter s suggesting and exploring binning as a solution. 
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FIG. 1: MG parameter evolution in redshift and scale modeled using a new hybrid method. We plot here a 3D representation 
for an example of the new hybrid binned evolution for the MG parameter Q(k,a) as given by our Eqs. (|12|) and (|13[) for the 
parameters Q(k,a) with Qi — 1.20, Q2 = 1.15, Q3 = 1.05, Q4 = 1.10, ztgr = 2, and k c = 0.01. We can see along the 
z-axis how the binned aspect can allow for different best fit values for the MG parameters in the redshift space while along 
the fc-axis we can see the monotonic evolution in k evolving from some large scale (small k) value to a small scale (large k) 
value exponentially. The hybrid parametrization combines the z-binning method that was shown to be robust with a smooth 
evolution in k space. 



To take advantage of all these techniques, we have developed two versions of our code. One version of the code 
uses the functional form (TTTT) . It provides the option to apply the functional form (|TT|) to either Q and R (as done in 
[5fjj ) or Q and T>. The other version of the code is based on binning methods. It provides the option of the second 
approach with traditional binning (described above), or alternatively the third, hybrid approach with two redshift 
bins, but the evolution in scale evolves monotonically. 

In the bin ning version of the code, we evolve only Q and V. Transitions between the redshift bin are evolved 
following [4J, [53[ and use a hyperbolic tangent function with a transition width z tw = 0.05. In this way the binning 
can actually be written functionally as (with X representing Q or T>) 

X (k, a) = 1 + Y k) + X ^ k) - X ^ k) tanh Z —^ + 1 ~ X »W tanh £Z£™ (12) 

Z Z Ztw ^ Ztw 

where Zdiv is the redshift where the transition between the two redshift bins occurs and ztgr is the redshift below 
which GR is to be tested. We hard code z TGR = 2z div to give us equally sized bins, but this of course is optional and 
can easily be changed. X Zi (k) represents the binning method for k in the ith z bin. For the suggested hybrid method 
it has the form 

X Z1 (k) = X ie - k l^ +X 2 (1- e- k ' k ") (13) 
X Z2 (k) = X 3e - fe ^+X 4 (l-e- fc /H 



while with traditional binning in principle evolves as 



X\ if k < k c 



X ^ = < x 2 iffc>C ( 14 ) 



Y (h\ - / Xs if k < kc 
Z2[ ' \X 4 if k > k c 



Here though, we have rather chosen to implement the traditional binning method with some control on the transition 
as: 

X 2 + Xi X 2 - Xx k - k c 
X Zl {k) = + tanh -^T — ( 15 ) 

v n\ Xx + X-i Xi-X-i k-k c 
X Z2 (k) = 1 tanh ■ 



ki 
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FIG. 2: We compare the effect the evolution of the MG parameters has on the matter power spectrum using both the binning 
and the new hybrid methods. The solid black line is the matter power spectrum produced using the best fit WMAP7 parameters. 
The blue dashed line was produced using the traditional binning method, and the red dash-dotted line was produced using the 
new hybrid method. The hybrid method produces a more realistic, smooth matter power spectrum and thus is physically well 
motivated. These plots were produced using parameters Qi = 1.20, Q2 = 0.88, Q3 = 1.52, Q4 = 0.95, Oi = 1.03, T>2 = 
0.95, X>3 = 0.99, T>4 — 0.93, ztgr = 2, and k c = 0.01. The amplitudes of the two modified matter power spectra are 
normalized to the GR amplitude at k ~ 7 x 10 -2 



where k tw is the transition width between k bins. We set k tw = k c /10 since we want to imitate traditional binning 
and ensure the transition between the bins is very rapid compared to the scale at which the transition occurs, thus the 
evolution maintains the appearance of being described by actual bins while having a functional definition. Although 
k tw can be used to smooth the transition between the bins but doing so to remove the appearance of actual bins also 
hinders the users ability to control the value of the MG parameters at k = 0, as is easily done in the hybrid method. 
We leave it to the users to choose their favorite method. 

To further illustrate the hybrid evolution, in Fig. Q] we plot an example of the evolution of the MG parameter 
Q(k, a) using this new hybrid form with Q\ = 1.20, Q2 = 1.15, Q3 = 1.05, Q4 = 1.10, ztgr = 2, and k c = 0.01. 
In this figure, one can see the binned aspect of the MG parameter evolution in redshift space, while the evolution in 
scale, fc-space, evolves monotonically. The advantages to this smooth evolution in scale are illustrated further in Fig. 
[3J where we compare the matter power spectrum produced using an identical set of parameters. The hybrid method 
produces a much more physical and smooth matter power spectrum. 

III. MODIFIED GRAVITY GROWTH EQUATIONS AND NUMERICAL IMPLEMENTATIONS 
A. Synchronous Gauge Variables and CMB Implementation 

We will first focus on our modifications to the publicly available code CAMB (Code for Anisotropics in the Microwave 
Background) [60( which calculates the various CMB anisotropy spectra (Cj T , Cf E , Cf E , Cf B ) as well as the three- 
dimensional matter power spectrum Pg(k, z). These observables are very powerful in constraining both the growth 
history of structure in the Universe as well as the expansion history of the Universe, and thus are very useful in 
constraining the parameters we use to test general relativity. 

CAMB is written in the synchronous gauge. Instead of using the metric potentials <fi and tp of the conformal Newtonian 
gauge, it uses the metric potentials h and 77 consistent with the notation of Ma and Bertschinger [65[. From Eq. (18) 
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of [65[ , we know the metric potentials in the two gauges are related to one another by 

<j> = rj — Ha, 
-0 = a + Ha, 
h 



where k a = 



37). 



(16) 
(17) 

(18) 



CAMB evolves the metric potential rj (actually evolving kr/) as well as the matter perturbations, 5i, heat fiux,g.;, and 
the shear stress <7j for each matter species in the synchronous gauge according to the evolution equations given in 
[65j . Additionally, the CAMB variables a cam b and Z are evaluated at each time step. They are useful for defining the 
evolution of the matter perturbations and are defined as 



ocamb = ka = 
h 



H 



2k- acAMB ~h- 



(19) 
(20) 



Using these variables allows CAMB to be written in such a way that the evolution of all other variables is changed 
simply by adjusting the evolution of the metric potential rj. Thus it is important that we derive an equation for the 
evolution of 77 consistent with the modified growth equations ([8]) and (fT7J|) . To do this, we begin by subbing (fT6)) into 
p]) and taking the time derivative. This gives 



Pi(a)(QAi + QAi) + QAi—pi(a) 

dr 



Ha + Ha 



where p~i{a) — 8irGa 2 pi. Next subbing (JSJl and (IT71) into (TT0|) gives an expression for a: 

a = -Ha - J2 Ma) [(22? - Q)A t + 3Q(1 + w^n] . 

i 

Now subbing Eqs. ([7]), as well as the time derivative of pi from matter conservation into (|2"TT) we get 



2k 2 



M^) [QAi - HQAi + 3Q(1 + wi) + Hip) - Qfij + 2HVA, - HQA t 



- (H 2 - H)c 



(21) 



(22) 



(23) 



where 

A =k 2 +3(H 2 -H). 

Next we can sub in for ip and <j> in Eq. (|2"3"j) by using (JTTJ) and the time derivative of (TTtl)) . This gives: 



-1 



f) = ^2 I] Ma) QAi - HQAi + 3Q(1 + w t ) (7) + (H 2 - H) 



(24) 



Qfx^+mVAi-HQAi -(H 2 -H)a. (25) 
k J 



We want an equation where all the variables on the right-hand-side are in the synchronous gauge, as those are the 
variables evolved by CAMB. We already know that A and a are gauge invariant, so we need not worry about them. 
However the qi in the above equations is still in the conformal Newtonian gauge. It transforms to the synchronous 
gauge according to Eq. 27b of 65] which by converting, as described above, from 9 to q can be written as 



4 C) =4 S) + (l + Wi)ka. 
So finally subbing ([26]) into (f25j) and collecting the terms we have [equivalent to Eq. A8 in [5C 



V= ^^2(H 2 -H)k 2 a + J2Ma) [( 



2H [2? - Q] + Q) At - Q(l + w t )k z a - Qf^ 



(26) 



(27) 



with 



f Q = k 2 + -Q^p l {l + w l ). 



(28) 
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After the metric potentials and matter perturbations are evolved the next major change is to redefine the deriva- 
tives of the Newtonian metric potentials, cb + ip, which go into evaluating the ISW effect in the CMB temperature 
anisotropy. We can get this quantity quickly by applying our knowledge that the quantities Aj and <Ji are invariant in 
transformations between the synchronous and conformal Newtonian gauges. Thus simply taking the time derivative 
of (fTO]) and subbing in for A and pi we get: 



{(l+3wi)Q + 2V)M-Q 



3(1 + w l )a i 3Q(1 + Wi)&i 



(29) 



-CDH - i>)Ai + 23(1 + wt) (k 2 a - 3r)) + 2?/i 



Other small changes to the scalar source term are necessary and were calculated using a slightly modified version of 
the Maple worksheet lineofsight. txt available from 60] . 

Some additional modifications though still remain to alleviate problems caused by the likely possibility that bad 
combinations of MG parameters will be introduced to CAMB by the sampling routine of CosmoMC. We encountered 
three routines which routinely ran into problems with bad combinations of MG parameters, the reionization module, 
the halofit module and the CMB lensing module. To alleviate these issues we added error variables to these routines 
where in the event a bad parameter combination causes these modules to fail, the calculation is aborted and that 
parameter combination is rejected rather than causing the CosmoMC to stop. With all of these changes in place CAMB 
is able to calculate and output all the modified CMB spectra as well as the modified matter power spectrum. 



B. Weak-Lensing Implementation 

Weak-lensing tomography shear-shear cross correlations are useful to constrain both the expansion history and the 
growth history of structure in the Universe. As discussed above this data can give us direct measurements of the 
parameter 23,;thus it is important to have a module for CosmoMC which takes advantage of this data. Here we discuss 
our module using the weak-len sing tomography shear-shear cross correlations of the HST-COSMOS survey recently 
compiled by Schrabback et. al. [6J|. It is worth noting here that this data set probes scales on which perturbations are 
not described by linear theory alone. The parametrized modified growth equations (JSj) and (|10[) describe deviations 
from GR on linear scales only. Though we use the halofit module to describe the nonlinear part of the power-spectrum, 
it should be noted that this may not be completely accurate and should be considered a source of possible systematic 
error for our MG parameter constraints. The authors of [63| perform a refined analysis of the HST-COSMOS survey 
of (66|, in combination with the COSMOS-30 photometric redshift catalog provided by 67]. The shear-shear cross 
correlations were calculated between 6 redshift bins, 0.0 < z < 0.6, 0.6 < z < 1.0, l.Q < z < 1.3, 1.3 < z < 2.0, 
2.0 < z < 4.0 and a sixth bin that contains all faint galaxies with a numerically estimated redshift distribution from 



0.0 < z < 5.0; see Fig. 6 of [63j. Certain exclusions were made as described in [63[, such as the luminous red galaxies 
(LRGs) to avoid G-I intrinsic alignment bias and the lowest angular theta bin because of model uncertainties. Only 
bright galaxies (i < 24) were used in the first bin (z < 0.6). Also, autocorrelations were not used in bins 1 — 5 to 
reduce the effect of I-I intrinsic alignments. We modified the publicly available code for the COSMOS 3D weak-lensing 
built by Lesgourgues et. al. 68] to incorporate the shear cross correlations as seen in, for example, [63| and the MG 
parameters. 

The shear cross correlation functions £v _ (8) between bins k, I are given by 

1 r°° 

= ^J dee j 0A (te)p* l (£), (so) 

where J n is the n*' l -order Bessel function of the first kind, £ is the modulus of the two-dimensional wave vector, and 
P^ 1 is the convergence cross-power spectra between bins k,l. 

In order to evaluate the convergence cross-power spectra, let us first recall that working with the metric ([IJ, the 
convergence is defined as 

K (n) = J** d X g{ X ) !k(x) \"J 2 [4>(f K ( X )n,x) + 1>(Jk(x)*, x)} , (31) 

with comoving radial distance X , comoving distance to the horizon comoving angular diameter distance Jk(x)^ 
and the weighted geometric lens- efficiency factor g( X ) given by 
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corresponding to the galaxy redshift distributions p. Applying the Limber approximation to (|31[) to get the convergence 
cross power spectrum in the absence of any modifications of growth gives us the usual result in terms of the three- 
dimensional (nonlinear) power spectrum Pg: 

where gk{x) is the weighted geometric lens-efficiency factor as defined above corresponding to the k th redshift bin. 

Now to incorporate our MG parametrization we substitute the late-time form (ie. negligible anisotropic stress) of 
(|f OP into (|31[) and again apply the Limber approximation. This gives us a modified convergence cross power spectrum 
that takes into account our MG parametrization: 



In calculating our cross correlation functions we follow instructions from [63| and weight the predictions for the k th 
bin, 6k, with logarithmic spaced upper and lower limits 9k,max and 9k, m in, respectively, based on the number of 
galaxy pairs, N for a given 9 as 

N(6')£ + -(6')d6' 

£(0k) = (35) 

N{e>)d6> 

with N{9) oc 0(0.0004664 + (9(0.0044118 - 8.90878 x lO" 5 !?)) 0, 

To evaluate the likelihood we use the 160 x 160 dimensional covariance matrix for the corresponding weak shear 
correlations provided by Schrabback et. al. [H[ applying the correction to the inverse covariance, C _1 , seen in |63ll70j. 

as 

C _1 = 0.4390 C -1 . (36) 

for 288 independent realizations and a 160 dimensional data vector. We then account for the 10% uncertainty in the 
numerical estimate of the galaxy redshift distribution of bin six, pq(z), by marginalizing over the nuisance parameter 
f z which is applied following [63] where for the sixth bin we instead use the galaxy redshift distribution p' e (z, f z ) given 
by p' e (z, f z ) = Pe(f z z). Now the likelihood (£) is given by: 

-21n£ = ^(^.-U s ) 4 [c- 1 ] ij (6/ 1 -U s ) J + (^ T ^) ■ (37) 



C. ISW-Galaxy Cross Correlations Implementation 

As discussed in, for example, [H, [H3| , cross correlations between the ISW effect and the galaxy density can be 
very useful in constraining the MG parameters. Thus, we modified the publicly available module by Ho et al. [6ll [62| 
for calculating these cross correlations in the presence of a modification to gravity described by Eqs. (|SJ) and (|10[) . 
Let us review some of the details of calculating these cross correlations. 

First recall, when working with the perturbed FLRW metric (UJ that the CMB temperature anisotropy due to the 
ISW effect can be written as [6l| 

AT ISW (h) =J °d T ^-^ + ^), (38) 

where r r is the conformal time of recombination and to is the conformal time today. The angular galaxy density 
fluctuations that we want to cross correlate these temperature fluctuations with are 



g(n) = / dz b(z) U(z) S m (f K (x(z))h, z) , (39) 



where b{z) is the galaxy bias, H(z) is the normalized selection function, and fK(x( z )) is t ne comoving angular diameter 
distance redshift z. 
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Following [6l|, we can once again use the Limber approximation to write the cross-power spectrum between these 
two fields in the absence of any modifications to the growth as 

C * = cHl+ l/2)i J dzb(zMz)—D(z)- [(1 + z)D(z)] P [j^^ ) , (40) 

where the growth factor D(z) = D(<S)8 m (k, z)/S m (k, 0) [please note that this D(z) is not the parameter V we use to 
test GR] is calculated using the standard differential equation for 8 m (z) at linear scales, S m + 2H8 m + k 2 ip/a 2 = 
with ip given in Eq. Pdj) . This result was attained by considering at late times in GR the potentials in (|38p can be 
written combining to ^ and © as 

4>(k,z) = i,(k,z) = -^(l + z)C2 m ^M. (41) 

When testing general relativity though, Eq. (f4"Tj) no longer holds and we must instead, as with the weak-lensing 
cross correlations, refer to the late-time version of (|T0|) giving 



<f>(k, z) + ip(k, z) = -iHl T>(k, z) (1 + z) n J mi ^ Z \ (42) 
With these changes the ISW-galaxy cross power spectrum is given by: 



Additionally, to fully account for the gravity modifications, the differential equation used to calculate the 8(k, z) in 
D(z) must be modified. Combining Eqs. © and ([TUf to get k 2 ip, the late-time form of differential equation for 8(k, z) 
now reads 

S m (k, z) + 2H5 m (k, z)-^- [2V(k, z) - Q(k, z)] (1 + z) 3 Q m S m (k, z) = 0. (44) 

In order to correctly perform these cross correlations, the code by Ho et al. also corrects factor b{z)H{z) to account 
for a magnification bias on SDSS quasars and NVSS radio sources [61j . In calculating this magnification bias, the 
weak-lensing convergence must be considered on scales smaller than the Hubble distance [7l|. Thus, we modify the 
code to account for changes in the weak-lensing convergence according to (|3Tj) and (|T0|) but consider only large-A: 
values of the MG parameters for these corrections in calculating the magnification bias. 

The original code by Ho et al. also has a likelihood code which uses cross correlations between weak-lensing of the 
CMB and the galaxy density. A detailed description of this part of the code is contained in [62j|. We modify this 
portion of the code to account for the changes to the convergence in the presence of modifications to the growth as 
described by Eqs. (HDJ and (|3T1) . 



D. WiggleZ Baryon Acoustic Oscillation Likelihood 

Recently, the WiggleZ Dark Energy Survey released its full data set of BAO measurements (64|. We include a 
likelihood module for this data set in ISiTGR. Let us quickly review the likelihood calculation for this data. 

The data set released by [64| includes measurements of the acoustic parameter A(z) at three different effective 
redshifts: z = 0.44, z — 0.6, and z — 0.73. The acoustic parameter was first introduced by [72j and is defined as 



= ioozMz)VCT (45) 

cz 



where O m /i 2 is the physical matter density. Dy(z) is the effective distance defined by: 

D v {z)={D\(z){l + zf^->j / \ (46) 

with the angular diameter distance to the redshift z, Da{z). 

To calculate the likelihood, the theoretical value of A(z) is calculated at each of the effective redshifts. Then using 
the covariance matrix for the data given in (63 |. the likelihood for a given model is expressed as 

- 21n£ = (A* ~ A obs ) l [C- 1 ]^ (A th - A obs ) j . (47) 
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FIG. 3: We plot the 68% and 95% C.L. constraints on the MG parameters Qo, Ro, and T>o from evolving the parameters 
with a functional form. We use all available data sets included in ISiTGR: SN, BAO, AGE, H , CMB, MPK, ISW, and WL. 
All constraints using this method are fully consistent with GR, however, as we show further, the correlations between these 
parameters and some cosmological parameters are significant. 
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FIG. 4: Left panel: We plot the 68% and 95% C.L. constraints on the MG parameters Qi and T>i, i = 1,2,3,4 from using 
traditional bins for k and z. z-bins are < z < 1 and 1 < z < 2 with GR assumed for z > 2 and fc-bins are k < 0.01 and 
k > 0.01. These constraints come from using all available data sets included in ISiTGR: SN, BAO, AGE, H , CMB, MPK, ISW, 
and WL. While all constraints are consistent with GR, bin 2 seems to be indicating some tensions with GR values just within 
the 95% C.L. Right panel: We plot the 68% and 95% C.L. constraints on the MG parameters Qi and T>i, i — 1,2,3,4 from 
using the new hybrid binning method, z-bins are still < z < 1 and 1 < z < 2 with GR assumed for z > 2, while the transition 
scale for k evolution is k c = 0.01. Again, we use all available data sets included in ISiTGR: SN, BAO, AGE, Ho, CMB, MPK, 
ISW, and WL. All constraints using this method are fully consistent with GR but as we show further, the correlations between 
these parameters and some cosmological parameters are significant. 



IV. CORRELATIONS OBTAINED BETWEEN COSMOLOGICAL PARAMETERS AND MODIFIED 

GRAVITY GROWTH PARAMETERS 

We provide here our results about various constraints and correlation coefficients found between core cosmological 
parameters and modified growth parameters. We use the standard definition for the correlation coefficient: 

Corrip^) = (48) 

where p Xl p y are the parameters, Cov(p x ,p y ) is the covariance of the two parameters, and (j(p x ) and <j{p y ) are their 
respective standard deviations. 

In addition to varying a given set of modified gravity parameters, we vary the six core cosmological parameters: 
Ct^h and the f2 c /i 2 , the baryon and cold-dark matter physical density parameters, respectively; 9, the ratio of the 
sound horizon to the angular diameter distance of the surface of last scattering; T re i, the reionization optical depth; 
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Correlation coefficients between U m and a$ and the MG parameters 


MG parameters evolved using traditional 


binning 


Parameter 


Qi 




Qz 


Q 4 




v 2 


Vz 


£>4 




0.0116 


-0.0535 


-0.0338 


-0.0371 


0.0114 


-0.0988 


0.0190 


-0.0818 


erg 


0.0937 


0.1540 


-0.1583 


-0.4267 


0.2093 


0.2857 


0.2377 


0.0185 




MG parameters evolved using hybrid binning 


Parameter 


Qi 


Q 2 


Qz 


Qi 


Th 


v 2 


Vz 


£>4 




0.1571 


-0.1639 


-0.003 


-0.0105 


0.1814 


-0.2172 


-0.0056 


-0.0070 


erg 


0.1692 


0.4282 


-0.1076 


0.1676 


0.2147 


0.5390 


0.1937 


0.0716 



MG parameters evolved using the functional form 



Parameter 



OS 



Qo 



-0.1173 



-0.5409 







-0.0294 



-0.2631 



Ro 



0.0542 



0.5742 



TABLE I: We list the correlation coefficients for O m and as and the various MG parameters. We see that the MG parameters 
for all evolution methods are significantly more correlated with as than Q m , although some correlation with Q m does exist. 
While traditional binning shows overall the least amount of correlation between the MG parameters and the two cosmological 
parameters, it, as shown in Fig. [2] has a problematic jump in the matter power spectrum. The hybrid binning is the next least 
correlated of the methods, with parameters from the functional form evolution showing the most amount of correlation with 
the two cosmological parameters compared here. 




FIG. 5: We plot here the 2D confidence contours for f2 m and erg and the MG parameters Qo and X\>, from using the functional 
form to evolve the MG parameters. As seen in Table [I] this evolution method overall has the most amount of correlations 
between the MG parameters and the two cosmological parameters. 



n s , the spectral index; and lnl0 10 A s , the amplitude of the primordial power spectrum. Additionally, the results 
illustrated here always use the following data sets for expansion history constraints, WiggleZ BAO measurements [13] , 
the supernovae Union2 compilation of the Supernovae Cosmology Project (SCP) [73| and references of other compiled 
supernovae (SN) therein. We also use the prior on Hq — 74.2 ± 3.6 km/s/Mpc given by [llj], and a prior on the age 
of the Universe (AGE) 10 Gyrs<AGE< 20 Gyrs. 

We include in our analysis the three parametrization approaches: functional form, binned method, and hybrid 
method. 

First, we utilize the functional form to evolve the MG parameters. For simplicity, we choose to assume the 
modifications are scale-independent (k c — oo) and evolve the parameters Q and R with the functional forms described 
by (HI]): 

Q(k,a) = (Qo-lK + 1, (49) 
R(k,a) = (,R -l)a s + l. (50) 

So, in addition to varying the six core cosmological parameters listed above, we vary the MG parameters Qo, T>q 
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(inferring i?o from Qo and T>q), and s. To constrain these parameters we use all of the available data sets: the WMAP7 
temperature and polarization spectra (CMB) [75j], the matter power spectrum (MPK) from the Sloan Digital Sky 
Survey (SDSSJDR7 [zH], the ISW-galaxy cross correlations [Hll E3], and the refined HST-COSMOS weak-lensing 
tomography [63]. 

Second, we bin the MG parameters rather than evolve with some functional form. We bin traditionally in both z 
and k as described by (fl"2"|) and (fl"5j) with Zdiv — 1-0, ztgr — 2.0, z tw = 0.05, k c — 0.01 and k tw = 0.001: 



Q(k, a) 



1 + 



. l+V z Ak) V z Jk)-V z Ak) z 
V(k,a) = 7 p^- + 2 w — ^^-tanh 



l_Q fi (fc) tanhi -2 

0.05 2 0.05 ' 

1 l-V z Jk) , z-2 



(51) 
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FIG. 6: We plot here the 2D confidence contours for Q. m and erg and the MG parameters Qi and T>i, i = 1, 2, 3,4 from using 
traditional bins for k and z. As seen in Table U this evolution method has the least amount of correlation between the MG 
parameters and the two cosmological parameters of the three evolution methods, but as discussed earlier suffers from the 
appearance of a jump in the matter power spectrum due to the rapid transition between scale bins. 
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FIG. 7: We plot here the 2D confidence contours for Q. m and erg and the MG parameters Qi and T>i, i = 1, 2, 3,4 from using 
the new hybrid method to evolve the MG parameters. As seen in Table [I] this evolution method has only a moderate amount 
of correlation between the MG parameters and the two cosmological parameters. 



with 



QzAk) 
QzM 



Q2 


+ 

2 


Qi 


h Q 2 


2 


Qi 


tanh 


k - 0.01 
0.001 


Qi 


+ 


Q 3 


Qi 




Q 3 


tanh 


fc — 0.01 




2 






2 




0.001 


Z>2 


+ 


Vx 


y v 2 




Vi 


tanh 


k - 0.01 




2 






2 




0.001 




+ 








v 3 


tanh 


A; -0.01 




2 






2 




0.001 



(52) 



(53) 



So, the z-bins are < z < 1 and 1 < z < 2 and GR is assumed for z > 2, while the fc-bins are k < 0.01 and k > 0.01. 
We vary all eight MG parameters, T>i ,Qi i = 1, 2, 3, 4, in addition to the six core cosmological parameters and again 
we use all of the available data sets: CMB, MPK, ISW, and WL. 

Finally, we evolve the MG parameters using the newly introduced hybrid method with evolution with redshift (z) 
dependence identical to that of traditional binning method, Eq. (|5"Tj) . and scale (k) dependence described by Eq. (Q2 
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with k c = 0.01: 



Q zi (k) 
V Z2 {k) 



Q ie - fc /°- 01 

Q 3 ^ fc/0 - 01 
X, ie -*/o.oi 

P 3e -fe/o.oi 



+ Q 2 (i 

+ g 4 (i 

+ 2? 2 (1 
+ 2? 4 (1 



e 



e 



e 



fc/0.01 



fc/0.01 



fc/0.01 



fc/0.01 



), 
), 



(55) 



(54) 



So, just as with traditional binning, the z-bins are < z < 1 and 1 < z < 2 and GR is assumed for z > 2 while the 
scale dependence transitions exponentially between two values for each redshift bin with a decay constant k c = 0.01. 
We again vary all eight MG parameters, T>i ,Qi i = 1, 2, 3, 4, in addition to the six core cosmological parameters and 
use all of the available data sets: CMB, MPK, ISW, and WL. 

We plot the 68% and 95% 2D contours for the MG parameters in Figs. [3] and 0] for the three methods respectively. 
We find that the constraints on modified gravity parameters are consistent with general relativity for all methods; 
however, the correlations between these parameters and some cosmological parameters are significant as we describe 
below. 

In TableU we report the results obtained for the correlation coefficients between the modified gravity parameters and 
cosmological parameters. The table shows that modified gravity parameters are significantly correlated with as and 
mildly correlated with Q m . For example, the degeneracies between as and modified gravity parameters are found to 
be substantial, especially for the functional form evolution method where Corr(Qo, as) = —0.54, Corr(Ro, as) = 0.57, 
and Corr(T>Q, as) = —0.26. For the hybrid method, about half of the bins seem to present some strong correlations, 
with the strongest being Corr(T>2, as) = 0.54. The traditional binning method seems to present less correlation than 
the hybrid method, but again about half of the bins have strong correlations with the largest being Corr(Q4, as) = 



Correlations between MG parameters and as have been discussed briefly by [52j. They too found that the correla- 
tions were significant. We find that for both the hybrid and traditional binning method that correlations between as 
and the T> parameters are positive meaning higher values of T> gives rise to higher values of as- We also find that, with 
the exception of the parameter P4, the high k parameters are more strongly correlated than the low k parameters, 
agreeing with the conclusions of [52j that these correlations are strongly dependent on the matter power spectrum, 
weak lensing, and ISW data sets. The corresponding plots are given in Figs. [5J HI and for the functional form 
method, traditional binning method, and hybrid method respectively. Looking at Fig. [6] we can see that there is not 
much difference between the slope of contours for as and T>2 or D4, while the correlations coefficients given in Table 
U are much larger for T>2 this can be explained by the bulges in each of those contours. For T>2 this bulge is toward 
the higher values of that parameter and at higher values of as thereby increasing the positive correlation, while for 
T>4 this bulge prefers lower values of said parameter, but still higher values of as presenting an anticorrelation that 
washes out the underling positive correlation of the rest of the parameter space. These bulges are caused by parameter 
constraints that are not completely Gaussian, as also seen in Fig. 2) Non-Gaussianity in the constraints also accounts 
for the odd shapes of the contours in Fig. [5j showing the correlations for the functional form MG parameters. 



Testing general relativity on cosmological scales is a topic of great interest in cosmology today. It allows an 
exploration of the origin of cosmic acceleration, as well as testing possible extensions to general relativity. Future 
high precision data will allow us to test gravity on cosmological scales to a very high level of precision. In this 
work we focused on the presence of correlations between modified gravity growth parameters and core cosmological 
parameters finding that some of them can be very strong. We provided the underlying modified growth equations 
and their evolution. We implemented the time and scale dependencies of the MG parameters using a functional form 
method, a binning method, as well as, a new hybrid method. The hybrid method provides a smooth evolution in 
scale and a binned redshift (time) dependence that was shown to be more robust than time functional forms. We 
used the most recent available cosmological data sets including BAO measurements from the WiggleZ Dark En ergy 
Survey [64j as well as the Two-Degree Field and DR7 Sloan Digital Sky surveys (7(1 [77}; priors on Ho from [74| : 
the Wilkinson Microwave Anisotropy Probe (WMAP) 7- year CMB temperature and polarization spectra (75|; the 
supernovae Union2 compilation, which includes the 557 type-la SCP [73[ and references of other compiled supernovae 
therein; the MPK from the SDSS DR7 0; the ISW-galaxy cross correlations using the 2MASS and SDSS LRG 
galaxy surveys [HI, [62j; and the recently refined HST-COSMOS weak-lensing tomography analysis in 63]. We also 
described our implementation of the modified growth equations into a numerical framework, ISiTGR, an integrated 
set of modified modules for the pac kages CosmoMC [59| and CAMB [§(| to be used to this end. ISiTGR modifies multiple 
publicly available codes: CAMB '6Cj, CosmoMC (59[, the ISW-galaxy cross correlation code [6l|, [62j, and the original 



-0.43. 



V. 
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COSMOS weak-lensing code [68j]. When using and referring to this package, the original codes should also be cited, 
as well as the above data sets when used in a particular analysis. Our analysis found that strong correlations between 
G 8 and modified gravity parameters are present, in addition to mild correlations with other cosmological parameters 
such as m . The degeneracies between a% and modified gravity parameters are substantial for the functional form 
and also for nearly half of the bins in the hybrid and binned methods. These degeneracies will need to be addressed 
when using future high precision data to perform these tests. 
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